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(-H , Abstract. In solving a system of ordinary differential equations by an interval method the 

Ph| approximate solution at any considered moment of time t represents a set (called interval) 

containing the exact solution at the moment t. The intervals determining the solution of 
a system are often expanded in the course of time irrespective of the method and step used. 
The phenomenon of interval expansion, called the Moore sweep effect, essentially decreases the 
04 I efficiency of interval methods. In the present work the notions of the interval and the Moore 

~^ ' effect are formalized and the Infinitesimal Moore Effect (IME) is studied for autonomous 

f~^ ■ systems on positively invariant convex compact. With IME the intervals expand along any 

r~^ , trajectory for any small step, and that means that when solving a system by a stepwise interval 

numerical method with any small step the interval expansion takes place for any initial data 
fj ' irrespective of the applied method. The local conditions of absence of IME in terms of Jacobi 

' ^y^ ■ matrices field of the system are obtained. The relation between the absence of IME and 

^^, simultaneous dissipativity of the Jacobi matrices is established, and some sufficient conditions 

of simultaneous dissipativity are obtained. (The family of linear operators is simultaneously 
dissipative, if there exists a norm relative to which all the operators are dissipative.) 



U ■ INTRODUCTION 



For solving systems of ordinary differential equations different classes of numerical 
methods with guaranteed error estimation including interval methods are used. In 
solving a system by an interval method the approximate solution at any considered 
moment of time t represents a set (called interval) containing the exact solution at 
the moment t. The detailed account of interval methods can be found in monographs 
by R.Moore [1] and S.A.Kalmykov, Yu.I.Shokin, Z.Kh.Yuldashev [2]. 

As a rule, all kinds of rectangular parallelepipeds with sides parallel to coordinate 
axes [1,2] are used as intervals, less frequently - ellipsoids [3], balls of fixed norm 
[4,5] etc. 

One of shortcomings of stepwise interval methods is the following. The intervals 
determining the solution of a system are often expanded in the course of time irre- 
spective of the method and step used. The simplest example of strong expansion 
of intervals during a short time, belonging to R.Moore, is given in [1]. The phe- 



effect was investigated only for some particular systems and particular intervals [1]. 

In the present work the notions of the interval and the Moore effect are formalized 
and the Moore effect is studied for autonomous systems on positively invariant 
convex compact. 

Formally, one can get rid of the interval expansion for any globally stable system 
(i.e. such a system, any solution of which is stable according to Lyapunov). To 
demonstrate that, let consider a smooth autonomous system: 

l^/W (1) 

on the positively invariant compact B C M". Construct a metric p on the set B, 
assuming for any x E B,y E B: 

p{x,y) = sup||a;(t) -y{t)\\, 

where x{t), y{t) are the solutions of the system (1) with the initial conditions 
a;(0) = X, 1/(0) = y. This metric is constricting for (1), i.e. for any pair x{t), y{t) 
of the solutions of (1) with the initial conditions in B 

p{x{t),y{t)) < p{x{s),y{s)) at t > s. 

The metric p is topologically equivalent to norm if and only if the system (1) is 
globally stable in B. If one considers as intervals all balls of the metric p, then in 
a definite sense the Moore effect is absent. That is, there is no interval expansion 
when constructing the exact interval solution with any step h > 0. The exact 
interval solution of X{t) is defined in the following way: X{0) = Xq, where Xq is 
the initial interval with the centre at the point a;(0) = Xq; X{{n+l)h) is the minimal 
interval with the centre at the point x{{n + l)h) containing ThX{nh), where Tt is 
the transformation of the phase flow of (1) during the time t > 0. Indeed, the radius 
X{{n + l)h) does not exceed the radius X{nh) at any n. 

If the system is not globally stable, then metric is not topologically equivalent to 
the norm. It means that small, in usual sense, intervals became large in the metric p. 
This circumstance makes one refuse from consideration of similar metrics. Moreover, 
if the system (1) is absolutely unstable (for example, a system with mixing), then 
there is no reasonable way to get rid of the Moore effect. 

The decsribed method of elimination the Moore effect for globally stable system 
is non-constructive. This can be demonstrated as follows: for constructing the 
constricting metric p one must know all exact solutions of the system (1). But then 
it is unreasonable to solve the system numerically. We must have constructively 
verifiable conditions of absence of the Moore effect and a way of construction of 
corresponding intervals. This is what we deal with in the present paper. The 
conditions of absence of the Moore effect are of local character and formulated 
in terms of Jacobi matrices of the system. Except that the causes of frequent 
appearance of the Moore effect will be pointed out. 

1 Interval spaces and the Moore effect 
1.1 Interval Spaces 



Definition 1. We call the family J of convex compacts in R" the interval space 
(and its elements - intervals) , if it satisfies the following conditions: 

a) JI is closed with respect to multiplication by non-negative scalars: 

if W^ e J, a > 0, then aW = {ax \ x e W} G J; 

b) J is closed with respect to intersection: 

if Wi e J, 1^2 e J, then WinW2e J; 

c) JI is closed according to HausdorfF (i.e. in the Hausdorff metric); 

d) if W^ G J, W ^ {0}, then G riW. 

Remind [6] that the Hausdorff metric on the set of all compacts in W^ is introduced 
as follows: 

Ph{x, y) = Taax{m.ax min ||a; — |/||, max min ||a; — i/||}, 

where x,y are the compacts in M", ||.|| is a fixed norm in M". All Hausdorff metrics 
in M" are equivalent. 

Further on by \im.Hi^oo Wi we denote the Hausdorff limit of the sequence {VFi}^ 
at i ^ oo. 

Give several examples of interval spaces. 

Example 1. J is the set of all convex compacts symmetric with respect to 0. It 
satisfies all the properties from a) to d). 

Example 2. JI is the set of all symmetric with respect to rectangular paral- 
lelepipeds (including non-singular), i.e. sets of the form 

{x = {xi,. . .,Xn) eW : \xk\ < ak {k = 1, . . . ,n)}, 

where Ofc > (/c = 1, . . . , n). It satisfies the properties a), b) and d). 

Let now {Wjj^j^ C J, MuiHi^ooWi = W with a^ being ak, corresponding to 
Wi. 

If pniWi, W) < e, then W C Wi + Pe, Wi C W + Pe, where P, = {x G R" : 
\xk\ < e {k = 1, . . . ,n)} (here a norm in the definition of the Hausdorff metric is 
the /°°-norm). Then for any x & W 

\xk\ < flfc +^ (A; = 1, . . . ,n) 
is true and for any x E Wi 

\xk\ <]mi^ooO'k^; 

lim a^^' < hmj^ooOfc , 
i.e. there exist the limits 

dk = lim a^ {k = 1, . . . ,n). 

If x G W^, then 
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j— >oo 



Xk\ < lim a^ {k = I, . . . ,n). (2) 



Obviously, x*^*^ = biX E Wi, i.e. there exists such a subsequence of {a;*^*^}^]^ that 
2;W £ yy.^ X = \imi^ooX(i). Hence W = {x E W^ : \xk\ < Ofc (/c = 1, . . . ,n)}, i.e. 
W E I and the property c) is also satisfied. 

In constructing interval methods of solving different problems it is, as a rule, the 
considered interval space that is made use of [1,2]. 

Example 3. Let ||.|| be a norm in M", ^ = {a; G M" : ||a;|| < r}, where r > 0. Let 

J={r| r >0}, 

i.e. J is the set of all closed balls (further on we omit the word "closed") of the 
norm. All the properties from a) to d) are satisfied. These interval spaces are used, 
for example, in [4,5]. 

Example 4. The construction of example 3 can be generalized as follows. Let 

1 1.11 be the finite set of norms in M", 



II' 



(fc) 



{xeM'' : \\x\\^<rk {k = l,...,m)} 



where r^ > {k = 1, . . . , m) and 



(k) 
rk ■ 
l<k<m 



Let J = {Wri,...,r„ '■ ^fe > {k = 1, . . . ,171)}. Obviously, J possess the properties 
a), b), and d). 

Note that the same element of JJ can be associated with different sets of {r^}. To 

demonstrate that, let m = 2, sup^-^odl^L/ll^lli) ~ ^- Then [ = Wi^c', where C is 
any number not less than C. Also, even if one of r^ is equal to 0, then 

Wr„...,r^ = {0}. 

To each compact W C M" can be juxtaposed the set 

K(W^)}r=i : rk{W)=rrmx\\x\\f^ (/c = 1, . . . ,m). 

liWel then w = Wr,n...nWr^. 

Let now the sequence {VFjj^i converge according to Hausdorff to the compact 
W, with Wi E J for all i. Similarly to example 2, from the inclusions 

WiCW + i''\ WcWi + i''^ 

satisfied for each e > for all i > io{e) derive the existence of the limits: 

ffc = lim rk{Wi) (/c = 1, . . . , m) 

and conclude that 

i.e. W E ^, and the property c) is satisfied. 

Example 5. Let Q be a compact convex body without symmetry centre (for 
instance, a triangle in R^), G intQ. Assume 

J = {aQ : a > 0}. 

J possesses the properties from a) to d). 

Remark 1. Example 5 can be generalized. For this purpose it is necessary to 
consider compact convex bodies Qi, . . . , Qm, the interior of each of them contains 0, 



1.2 Dissipative Operators 

In this section the properties of the operators dissipative with respect to compact 
are studied. First, let remind some notations. 

The affine envelope of the convex set W is denoted by AffW, the relative interior 
W (the interior of W in AffW) is denoted by riW, the relative boundary of W (the 
boundary of W in AffW) is denoted by rdW. For the boundary of the set X we use 
the notation dX, intX - for the interior of X, coX - for the convex envelope of X. 
By the sum of the sets of X and Y from M" we mean the set {x + y : x E X,y E Y}, 
by I - the unit operator. 

Let introduce a new notion. 

Definition 2. The linear operator A in the space M" is called dissipative with 
respect to the family of sets {Wi,} C M" if every set Wu is positively invariant with 
respect to the system 

dx , 

- = A.. (3) 

In other words, every W^, is invariant with respect to the semi-group of the 
operators exp{At) (t > 0). 

Below we consider operators dissipative with respect to families of convex com- 
pacts. In particular, the operator is dissipative with respect to the families of all 
balls of some norm (for this, dissipativity with respect to only one ball is sufficient) 
if and only if || exp(At)|| < 1 at allt > 0. Thus, in this case we come to the known 
definition of dissipativity with respect to the norm [7] . 

The set of all operators dissipative with respect to {W^} is denoted by K{{W^}). 

Remark 2. If an operator is dissipative with respect to the family of compacts 
and the interior of at least one of them is not empty, then it is dissipative with 
respect to some norm. 

Indeed, any symmetric with respect to compact convex body is a ball of some 
norm (see, for example, [7]). Choose as a ball the following set: 

S = co{WU{-W)} (4) 

where W is any set of the considered family of {VTz^}, for which intW ^ 0. 

However, if W^ is a compact and the operator is dissipative with respect to the 
norm whose ball is S (4), then it does not yet mean that the operator is dissipative 
with respect to W (see also example 8). 

Remark 3. From the invariance of a family of compacts with respect to the 
linear operator follows the invariance of the Hausdorff closure (i.e. closure in the 
Hausdorff metric) of this family. Therefore from dissipativity of the operator with 
respect to the family of compacts follows the dissipativity with respect to Hausdorff 
closure of this family. 

Let W^ be a convex compact in M" with G riW. In this case AffW is a linear 
subspace, and if the operator A is dissipative with respect to W, then AffW is 
invariant with respect to A. Introduce the following functional on the subspace 
L{W) of the space L(R") (of linear operators in M"), consisting of the operators, 
with respect to which AffW is invariant: 

HwiA) = sup fiwiAx). (5) 



It is easy to see that A G K{W) if and only if 

fiw{exp{At)) < 1 

for all t > 0. 

In particular, the operator A G L{W) is strongly dissipative with respect to convex 
compactW if exists such e > that /ivK(exp(74t)) < exp(— et) at all t > 0. In general, 
the operator A is strongly dissipative with respect to convex compact W if and only 
li A + el e K{W) for some e > 0. 

If H^ is a ball of the norm ||.||, then strong dissipativity with respect to W means 
the existence of such e > that || exp(At)|| < exp(— et) for all t > 0. We come to 
the definition of stable dissipativity with respect to the norm [11, 12, 19]. 

Introduce in L(M") the following functional: 

( A\ V ^^w{I + hA)-l 
lw{.A) = hm 

In the case, when W^ is a ball of some norm (i.e. ^w is a norm), arrive at the known 
definition of the logarithmic Lozinsky norm [9, 10]. 

Lemma 1. The operator A G L(W^) is dissipative (strongly dissipative) with 
respect to W , if and only if the inequality 7h^(A) < (jwi.^) < 0) is satisfied. 

Proof. Sufficiency. The following inequality is obtained in [9] 

||exp(At)|| <exp(7(A)t) 

where 7(A) is the Lozinsky norm of the operator A, corresponding to the norm ||.||. 
By literal repetition of the reasonings from [9] (with a substitution of the norm by 
Minkovski functional), one can obtain the inequaliny 

/iH/(exp(At)) < exp(7H/(A)t) 

for allt > 0, from which immediatelly follows the sufficiency. 

Necessity. Evidently, 

/iH^(exp(At)) = ^w{I + At) + o{t) (t ^ 0). 
Therefore, 

7^(A) = ^hm . 

Let £ > 0. If /ivi/(exp(At)) < exp(— et) at allt > 0, then 

... , exp(— e/i) — 1 

which proves the necessity. The lemma is proved. 

Assign a relatively open convex cone Qx{W) to every point x G rdW according 
to the rule: y G QxiW) if and only if there exists such e > that 

X + ey E riW. 

Lemma 2. For strong dissipativity of A with respect to convex compact W it is 
necessarv and sufficient that for everv Doint x G rdW the inclusion 



be true. For dissipativity of A with respect to W it is necessary and sufficient that 
for every point of X G rdW the inclusion 



Ax e Q^{W) 

be true. 

Proof. Note that the operator A is strongly dissipative with respect to W if and 
only if there exists such to > that fj^wi^ + ^^o) < 1- Indeed, the existence of such 
to for a strongly dissipative operator follows immediately from the negativeness 
of 7ty(A). Conversely, if fiw{I + ^^o) < 1) then there exists such e > that 
Hwil + {A + e/)to) < 1. But then 7iy(^ + £l) < 0, the operator {A + el) is 
dissipative. It means that A is strongly dissipative. 

If the operator A is strongly dissipative with respect to W, then, according to 
the above, for each x G rdW there exists such t^ > that (/ + txA)x G riW. It 
means that the vector Ax belongs to to the cone Qx{W). 

Conversely, let the latter condition be satisfied. According to the hypothesis of 
the theorem and convexity of W, for each x G rdW there exists the only positive 
number s = s{x) such that {I+sA)x G rdW. Show that Sq = inix^rdw s{x) > 0. Let 
it be not so. Then there exists such a subsequence {xn}n=i that lim„^oo s(a^n) = 0. 
Choose from {xn} a converging subsequence {a;„'}. Let x = lim^^ooa^n'- For every 
n E N and for every e > 

[I + {s{xn')+e)A]xJ^W. 

Passing to the limit, obtain 

(J + eA)x ^ riW 

which contradicts the hypothesis of the theorem. 

Thus, So > 0. For any to G (0; Sq) is true iiy/i.! + ^^o) < 1; i-e- the operator A is 
strongly dissipative. 

If A G K{W), then for any e > we have AX — ex E Qx{W) (for any x G rdW), 
i.e. Ax G Qx- Conversely, if Ax G Qx-, then Ax — ex E Qx at any £ > 0, and A 
represents a limit point of the family of dissipative operators, i.e. A G K{W). The 
lemma is proved. 

Remark 4. Immediately from the Krein-Milman theorem [8] follows that it is 
sufficient to require from the lemma conditions that inclusions be satisfied not for 
all points x G rdW, but for extremal points of W only. In particular, if W is 
a polyhedron, then it is sufficient to test its vertices only. Thus, to elucidate the 
question about dissipativity (strong dissipativity) of the operator with respect to the 
polyhedron, one should test only the fulfilment of finite number of linear inequalities. 

Remark 5. In the proof lemma 2 we have used the obvious fact: the closure of 
the set K{W). 

One more fact follows directly from lemma 2. 

Lemma 3. The set K{W) is a closed convex cone. The cone of all strongly dissi- 
pative with respect to W operators coincides with riK(W) and with lJe>o(-^(^) ~ 
el). If {Wu} is a family of convex compacts with G riW^ for all u, then K{{Wu}) 
is a closed convex cone. 

Remark 6. If intW = 0, then intK{W) = 0. Indeed, if AffW is invariant with 
respect to the operator A,, then A + eA, 4 KiW) at e ^ 0. If intW ^ 0, then 



Definition 4. The operator A G K{W) is called stable (or roughly) dissipative 
with respect to W, if A G intK(W). 

Definition 4 generalize the definition of the stable dissipativity with respect to 
the norm [11, 19]. 

Pass to the consideration of operators dissipative with respect to interval spaces. 
Let find out for which interval spaces J the interior of the cone K{J) is not empty. 

Let l^ be a set of all compact convex bodies in M„. Fix some norm ||.|| in M" and 
assume 

d(W) = min ||a;||. 

x&dW 

Lemma 4. The function d{W) is continuous according to Hausdorff on the set 
V. 

Proof. First note that ii X e V, Y eV, then pnidX, dY) < ph{X, Y). Indeed, 
let ph{X, Y) < e. Then X C Y + S^ where S'^ = {x G M" : ||x|| < e}. Let, further 
on, there exists such i/q G (dY) fl X that yo ^ S^ + dX. Construct at the point 
I/O a tangent hyperplane L to Y. Let / be the direction of the external normal to 
dY at the point yo orthogonal to L. Draw a ray from the point yo in the direction 
of / to the point xq of crossing with dX. Construct such a ball S of the norm {|.|| 
with the centre at the point x that yo G dS. The radius of S is larger than e and 
S nY = {yo}- Thus, if one constructs a ball S" C S* of the radius e with the centre 
at xo, then 

S'nY = 0. 

But then xq ^Y + S^, i.e. X (^ Y + S^ what is contrary to the assumption. 

The existence of such yo C dY that yo ^ X U {dX + S^) is also impossible, since 
then yo ^ X + S'^, i.e. Y (t X + Se- Consequently, dY C S^ + dX, and that means 
d{X) < d{Y) + e. Similarly, d{Y) < d{X) + e. It means that \d{X) - d{Y)\ < e, 
and the function d{W) is continuous on V. The lemma is proved. 

Lemma 5. For non-emptiness of intK{J) it is necessary and sufficient for all the 
elements of the interval space J, exept {0}, to posess non-empty interior. 

Proof. Necessity. Follows immediately from remark 6. 

Sufficiency. Show that under the conditions of the theorem the inclusion 

-/ G intK{J) (6) 

takes place. 

To each point x {\\x\\ = 1) we assign the set W^ according to the rule: 

W^= fl W. 

WBx,wes 

According to the conditions b) and c) from definition 1, Wx G J. The set 



W= [j Wx 

\\x\\=l 

is compact. Indeed, W is contained in any element of J containing unit ball of the 
norm ||.||; such an element exists due to non-emptiness of the interior of all intervals 
(exept {0}) and the property a) from definition 1. Note that Hausdorff closure of the 
family {W^ : ||x|| = 1} represents a compact in the Hausdorff metric, contained in 



of lemma 4) the existence of such e > that d{Wx) > e for all such x that ||a;|| = 1 
(indeed, d{Wx) > 0, since e intWx)- 

Thus, there exists such e > that for all x (||a;|| = 1) the inclusion 

Ax e intWx 

is true if \\A\\ < e. 

In other words. Ax — x (z QxiWx) if ||^|| < ^i II^^H = 1 (see lemma 1). The more 
so, as Ax — X E QxiW) for all W^ G J {W 3 x, \\A\\ < e) at all such x that ||a;|| = 1. 
But then Ax - x e Qax{aW) for all a > 0, ||A|| < e. Hence, A- I e K{I), i.e. (6) 
is satisfied. The lemma is proved. 

Thus, we have shown that under the conditions of lemma 5 K{J) is a convex solid 
cone. 

Definition 5. The operator is stable dissipative with respect to the interval space 
J if it belongs to intK{S). 

For stable dissipative operators the remark 2 is true: if an operator is stable 
dissipative with respect to the family of compacts and the interior of at least one of 
them is not empty, then it is stable dissipative with respect to some norm. 

1.3 The Moore Effect for Autonomous Systems 

The results of the previous section can be applied to the study of the Moore sweep 
effect. First give the exact definition of what we understand by the Moore effect. 

Let in the vicinity of a compact convex body 5 C M" be given a smooth au- 
tonomous system 

f^/W (7) 

with B positively invariant with respect to (7), and let x{0) be determined inexactly, 
namely 

x(0) exo + Wo, 

where xq G B, Wq G J, xq + Wq E B, J is some interval space (see definition 1). 

Remark 7. Irrespective of particular numerical method (i.e. dealing with the ex- 
act solution of the initial value problem for (7) with the initial conditions a;(0) = xq) 
a stepwise interval solution with step h > can be described as follows. 

Let Tfi be the transformation of the phase flow of (7) during the time t {shift over 
time t), Wq G J is the initial interval (its sense is an uncertainty in initial data). 
Assume 

Xo = xo + Wo, 

Xm+l = T(^m+l)hXo + Wm+1, 

W^+i = {^ w, 

WZ}Wm+x{h),W(^^ 
Wm+lih) = Th{TmhXo + Wra) - T(^m+l)hXo. 

The sequence {Xm}^o is the exact stepwise interval solution of (7). 

Definition 6. The absence of infinitesimal Moore effect (IME) means that 
for any h > the sequence {W^m}m=o ^^ enclosed: W^ D Wm+i ^or all m, i.e. the 
obtained intervals do not expand. 



any small step the interval expansion takes place for any initial data irrespective of 
the applied method (since it is true even for exact solutions). 

Generalizing the construction [10] for norms, introduce the following functional: 

,, / X ,. fJ'wix + hy) - /iH^(x) 

Literally (with substitution of the norm for Minkovski functional) repeating the 
reasonings from [10] (pp.127, 426), come to the following statemets. 

Statement 1. If x{t) with values in M" is differentiable on connected subset T of 
the real axis, and W^ is a convex compact (0 G riW), then the function fiw{x{t)) is 
almost everywhere differentiable on T and the derivative (where it exists) coincides 
with the right-hand derivative, equal to Nw{x(t),s(t)). The right-hand derivative 
of fiw{x{t)) exists everywhere on T except the right-hand end. 

Statement 2. 

7ty(A) = sup Nivix,AX). 
xew 

By f'{x) further on we denote the mapping derivative of /. 

The main part of further results on IME can be obtained from the following 
theorem. 

Theorem 1. Let in the region U C i?" be given a smooth autonomous system 
{7), B G U he positively invariant with respect to (7) compact convex body. IME 
is absent for compact B, system (7) and interval space J if and only if 

fix) e K{I) (8) 

for all X E B, i.e. for any x E B the Jacoby matrix of system (7) in the point x is 
strongly dissipative with respect to JJ. 

Proof. Sufficiency. Let W E^. Consider two solutions xi{t),X2{t) of system (7) 
with initial conditions from B. Denote A(t) = xi{t) —X2{t). Using statements 1 and 
2 and the theorem on finite increment, estimate the derivative of /ipi/(A(t)): 

-/iH'(A(t)) = iVp^(A(t),rfA(t)M) < 
< sup iVp^(A(t),/'(a;e(t))A(t))< 

O<0<1 

< sup -iwU'Mt)))iiw{^{t)), 

O<0<1 

where Xc{t) = Xi(t) + Q{x2{t)—Xi{t)), < < 1 for all t > 0. By (8) and statement 
1 we obtain 

j^l^wiMt)) < 0. 

Since the latter inequality holds for allt > and for all W E S, in systen (7) on 
B IME with respect to JI is absent. 

Necessity. Let W^ G JI, xq E intB, to > 0, y E Af fB, y j^ 0. There exists such 
ho > that xq + hoy E B. Due to smoothness of system (7) there exist and are 
unique the solutions Xi{t),X2{t) of the initial value problem for (7) with the initial 
conditions xi(to) = xq, X2(to) = Xq + hoy. Assume A(t) = xi(t) — X2{t). Then 



A(to) „, , A(to) 



Nw\ ^Z,„ f'ix,) 






where Xc = xq + Qhoy, < < 1. 
By virtue of absence of IME 

-\nfiw{A{t))<0 

for allt > 0. Since if Xq + hoy G b, then: 

(a) Xq + hy E B for all h G [0; ho], 

(b) a set of those h G [0; ho], for which 

\fj.w{y> fj^wiyjj 

is dense on the segment [0; ho], 

and (c) due to its closureness coincides with this segment. 

By virtue of arbitrarity of the choice of Xq for any Xo G intB, t > 0, y G AffB, 
y j^ the inequality 

\fiw{y) i^w[y)) 

is satisfied. It holds also for any x E B,t > 0,y E AffB, y ^ 0. Hence, from lemma 
1 and statement 2 immediately follows dissipativity of f'{x) with respect to J for all 
X E B. The theorem is proved. 

Definition 7. The family of linear operators {A^} is called simultaneously dis- 
sipative, if there exists a norm relative to which all the operators are dissipative. 

Simultaneously dissipative operators were studied in detail in [11, 12, 17-22]. 

From theorem 1, example 3, and remark 2 we obtain the following theorem. 

Theorem 2. For existence of interval space in which at least one interval posesses 
non-empty interior and with respect to which in system (7) there is no IME on B, 
it is necessary and sufficient for the family {f'{x) : x E B} to be simultaneously 
dissipative. 

Thus, the problem of existence of the interval space, with respect to which IME 
is absent, is reduced to the problem of simultaneous dissipativity of Jacobi matrices. 
As sought for space one can choose a set of all balls of that norm relative to which 
all Jacobi matrices are dissipative. This norm is constricting for (7) on B (i.e. the 
distance between two solutions with initial conditions from B will not expand with 
time). Hence, all systems without IME (with respect to some interval space) on B 
are globally stable in B (see introduction). 

Bellow by C^{B) we denote the Banach space of smooth mappings of B in MJ^ 
with the norm 

df 



^i(,)=max||/(a:)|| + ^max 

fc=i 



dxi 



where 11. 1 1 is a fixed norm in 



Immediately from lemma 3 and theorem 1 the following statement can be ob- 
tained. 

Theorem 3. The set of systems on B without IME with respect to J is closed 
convex cone in C^{B). 

For this cone we use the notation Fb(J). 

Further on, speaking about the vicinity of an autonomous system in C^{B) we 
mean a part of the vicinity, consisting only of those systems for which the set B is 
positively invariant. 

Let us study under what conditions the interior of the cone -Fb(J) is non-empty. 

Theorem 4. For non-emptiness of intFB{T} in C^{B) it is necessary and suffi- 
cient for all elements of JJ, exept {0}, to possess non-empty interior. 

Proof. Necessity. Let exist such a set W^ G J that intW = 0. Consider any 
system (7) without IME with respect to J on B. Since intB ^ 0, there exist 
two different concentrical balls 5*1 and 5*2 of usual /^-norm, belonging to intB with 
5*1 C 5*2. Construct such a function g G Coo(IR") that g{x) = 1 for all x E Si and 
g{x) = at all a; ^ 5*2. Since AffW ^ M", one can construct a linear operator 
A G L(M") mapping AffW into such a subspace Eo ^ {0} that {AffW)nEo = {0}. 

Consider the system 

= f{x)+eg{x)AX, (9) 

ay 

where e > is arbitrary. The set B is positively invariant with respect to (9), since 

the vector field generating (9) coincides with / in the vicinity of dB. On the other 

hand, there exist Jacobi matrices (9) relative to which AffW is not invariant, i.e. 

in (9) exist IME with respect to J on B. Since in any vicinity of / there is at least 

one vector field, generating (9), then 

intFBil) = 0. 

Sufficiency. Consider the system dx/dt = —x. It is a system on B without IME 
with respect to J. Furthermore, if all elements of J, except {0}, possess non-empty 
interior, then by lemma 5 the matrix of the system is stable dissipative with respect 
to J (see definition 5). 

Consider the system: 

JL = -x + v{x), (10) 

dy 

where llf IL tm < £ with e chosen so that 

A-I e K{S) 

if \\A\\ < e (see the proof of lemma 5). Then all Jacobi matrices (10) are dissipative 
with respect to J, and if v is choosen so that B is positively invariant with respect 
to (10), then in (10) IME is absent (by theorem 1). The theorem is proved. 

Thus, intFeil) 7^ if and only if intK^) ^ 0. 

It is easy to see that in the proof of sufficiency in theorem 4 one can instead 
of the system dx/dy = —x consider any system whose Jacobi matrices are stable 
dissipative with respect to J. 

Remark 8. By analogy with the space C^{B) one can construct the Banach 
spaces C^{B) {k G N) with the norm 



where a = (ai, . . . , «„) is a multiindex: 

dHf 



|q;| = tti + . . . + an,D"f 



dx"^ . . . dx^" 

and the metric space C°°{B) with the system of seminorms 

{max||(D„/)(x)||: |a| < m}+- . 

Small C'^-additions {1 < k < +00) are small and in the C^-norm. Therefore, for 
C'^-smooth systems under the conditions of theorem 4 the interior of -Fb(J) is non- 
empty and in C^{B). As shows the proof of theorem 4 (necessity), if conditions of 
the theorem are not satisfied, then the interior of -Fb(J) in C^{B) is empty. 

Let clarify what autonomous system without IME in specific interval spaces looks 
like. 

Theorem 5. Any system without IME with respect to J from example 1 has the 

form: 

dx 

—- = ax + c, 

at 

where a < 0, C E M" is a constant vector. 

Proof. Let A G K{J). All the segments symmetrical with respect to belong to 

J. Every such a segment has the form {y G MJ^\y = ax, \a\ < 1} for some a; G M". 

The cone Qx (see lemma 2) for each segment consists of vectors of the form ax, 

where a < 0. Thus, every non-zero vector a; G M" is eigenvector of the operator A, 

corresponding to non-positive eigenvalue. Thus: 

K{I) = {al\ a < 0}. (11) 

Let now a system without IME have the form 



dxi 
"dt 

HiJb Yl 

"df 
According to (11) and theorem 1 



fi{xi,...,Xn); 

Jn\Xi, . . . , Xfi) ■ 



= o(2^j); (12) 

< 0. (13) 



dxj 
dfi _dh _ _ df, 



dxi dx2 dxn 

From (12) follows that fk depends only on Xk {k = 1, . . . ,n). It means that 
dfk/dxk also depends only on Xk, i.e. by virtue of (13) dfk/dxk = const {k = 
1, . . . ,n). Then 

dfi _ -9/2 _ _ dfn 



dxi dx2 dxr, 

and the system has the form: 



= a<Q 



dx 



a.x 



Thus, whatever nonhnear (or even hnear with non-scalar matrix) system we con- 
sider, if we take as J the interval space of example 1 (or any wider space), IME will 
be present in the system. From theorem 5 also follows that any dissipative with 
respect to all norms operator has the form a/, where a < (see also remark 3). 

Example 6. Consider J from example 2. J contains all symmetrical with respect 
to segments of coordinate axes (thus, the conditions of theorem 4 are not satisfied, 
i.e. iniFsiS) = 0). Let A G K{S). Reasoning like in proof of theorem 5, conclude 
that all coordinate axes are eigenspaces of the operator A, corresponding to non- 
positive eigenvalues. In other words, the matrix of the operator A is diagonal and 
non-positive. On the other hand, by virtue of lemma 2, all such operators belong to 
i^(J). Thus, systems without IME with respect to J on i? have the form 

dxi 

- at 
where 

|^<0(A; = l,...,n) 

for all X E B. 

From the considered example follows that when using standard intervals (rectan- 
gular parallelepipeds) IME will be observed in almost all systems in M" if n 7^ 1. 

The systems without IME with respect to J from example 3 on 5 represent all 
systems for which the norm {|.{| is constricting in B (see the text after theorem 2). 

Remark 9. Note that testing of dissipativity (stable dissipativity) of the operator 
with respect to the norm is equivalent to non-positiveness (negativeness) of the 
corresponding Lozinsky norm. For some norms an explicit form of corresponding 
Lozinsky norm is known (see, for example, [9] or [10, p.463-465]). In particular, 
for the Euclidean norm the Lozinsky norm of the operator A coincides with the 
largest eigenvalue of the operator {A* + A)/2. The Lozinsky norm of the operator 
A represented by the matrix (oij)" -^^^ with respect to l"^— and /°°-norms is given by 
the formulae, respectively: 

max (Re an + y^ iQjil); 
max (Re an + 5" laj,]). 

l<i<n ~^ 

In remark 9 it is assumed that the operator A acts in the space C". The definitions 
and used here properties of dissipative operators in complex spaces are analogous 
to those in real ones. 

Example 7. Let JI be the interval space from example 4. Then 

l<fc<m 

where i^iMi is the cone of all operators dissipative with respect to the norm i^iMi . 



sets of an autonomous system. Similarly, 



mtir(J) = Pi intK II 

l<k<m 



Let, for example, ||.||i be /°°-norm, ||.||2 be /^-norm in M^. Then the conditions 
of stable dissipativity of the operator A with the matrix (ajj)^^^^ with respect to J 
according to remark 9 are of the form: 

4aiia22 > (ai2 + 021)^; 

On + |ai2| < 0; 

|a2l| + 1^22 < 0. 



Thus, for the system of the form 

dxi 

dx2 



if the inequalities 



4 



dhdh 
dxi 



dxi 
5/2 



+ 



0X2 



fi{xi,X2); 

f2{Xi,X2) 



^ (dj\^df2 



(14) 



dx2 dxi 



dxi 



+ 



dx2 
dh 

dX2 



< 0; 



< 



are satisfied and the compact convex body B is positively invariant with respect 
to (14), then in (14) IME with respect to J (from example 4) is absent on B. For 
example, such is the following system: 



—rr = -2xi + X2] 




= 2a;i — 3x2 

if i? is the square {(a;i,a;2) : |a;i| < 1, |a;2| < 1} or the circle {(a;i,a;2) : xf = X2 < !}■ 
Example 8. Consider J from example 5. Let Q be rectangular triangle with 

vertices at the points (—1; 2); (—1; —1); (1; —1). 

From lemma 2 and theorem 1 follows that the cone -Fb(J) consists of the systems 

of the form (14), with respect to which the compact B is positively invariant and 



for which 



dxi dx2 
dxi dx2 



dxi 



df2 



dxo 



-3^ + 6^-2^ + 4^^' 



dxi dx2 



dxi 



dxo 



3^ - 3^ + 2^ - 2^^' 



< 0; 

< 0; 

< 0; 

< 0; 

< 0. 



(15) 



dxi dx2 dxi dx2 
is true. 

Substituting all the inequality signs in (15) by strict ones, obtain intFB{I). For 
example, the system 

dxi 



dt 
dx2 
"dt 



—Qxi — 4^2 



belongs to intFB{I) for B = Q. 

Corresponding ball S (see remark 2) is the parallelogram with the vertices in the 
points (—1; 2); (—1; —1); (1; —2); (1; 1). From remark 2 follows that 

ir(J) = K{Q) c K{S). 

One can see that i^(J) 7^ K{S). For example, the operator given by the matrix 

-4 -1 

2 

is dissipative with respect to S, but it is not dissipative with respect to J. In other 
words, in the systems without IME with respect to {aS : a > 0} (i.e. constricting 
according to the norm whose ball is J) there can be observed IME with respect to 
J. 

This example can be generalized as follows. Consider J from remark 2. In system 
(7) on B IME is absent with respect to J if the operators f'{x) for all a; G -B are 
dissipative with respect to all sets Qk{k = 1, . . . , m). 

To sum up, one can say the following. When using sufficiently wide interval 
spaces in almost all systems in accordance with theorem 4, IME is observed. In 
particular, IME takes place almost for all systems when using standard intervals 
(see example 6). Expansion of the interval space results in the appearance of new 
systems with IME: thus, in using a set of all symmetrical to convex compacts 
IME is absent only for linear systems with non-positive scalar matrices. And the 
most impotent: the question about the existence of interval space, with respect to 
which in the considered system IME is absent, is reduced to the problem of joint 
dissipativity of the Jacobi matrices. Therefore, there is no interval space with respect 
to which all for even if in some sense almost all) eloballv stable svstems would have 



corresponding interval spaces for each particular system. These problems are solved 
constructively very rarely. 

We have treated the Moore effect in a very strong sense. The condition of bound- 
edness of the sequence of intervals {VFm}^'^Q at any step h > (see remark 7) 
is weaker (and acceptable, generally speaking, for constructing sufficiently narrow 
interval solutions). This condition can be called the condition of absence of the 
asymptotic Moore effect (AME). It is the weakest from acceptable conditions, since 
with AME it is impossible to use stepwise interval methods to obtain narrow interval 
solutions at large times. The study of AME is still not completed. It is evident only 
that for a linear autonomous system in considering the interval space from example 
3 AME is equivalent to IME. One can suggest a hypotesis: the problem of existence 
and constructing of the interval space with respect to which AME is absent in the 
autonomous system is reduced to the question of simultaneous dissipativity of Jacobi 
matrices (and of constructing a constricting norm). 

2 Conditions of Simultaneous Dissipativity of Operators 

2.1 Some General Results 

In the present section some conditions of simultaneous dissipativity of the operators 
will be considered (see definition 7). 

A definition of a simultaneous dissipativity can be generalized in such a way. 

Definition 7'. A family of linear operators {Aa} is called simultaneously stable 
dissipative if there exists a norm with respect to which all operators Aa are stable 
dissipative. 

Lemma 6. Let the space R"^ be expanded into direct sum of subspaces Ei {i = 
1, . . . ,k) and each of them is invariant with respect to all operators of the family 
{Aq,}. Further on, let restriction of the family {A^} on any E^ be simultaneously 
(sumultaneously stable) dissipative. Then {A^} is simultaneously (simultaneously 
stable) dissipative. 

Proof. Let ||.||^ {i = 1, . . . , k) be the norms in Ei in which the restrictions of 
{Aa} on Ei are simultaneously (simultaneously stable) dissipative. Define the norm 
in M" in this way: 

k 

II II V^ II II 
X = > rT- 

11-^11 Z^ l|-^«IU5 

i=l 

where x = Y.i=i ^i with Xi E Ei {i = 1, . . . ,k). 

In this norm all operators A^ are simultaneously (simultaneously stable) dissipa- 
tive. The lemma is proved. 

It is known [7] that for one operator the norm with respect to which it is dissipa- 
tive exists if and only if the spectrum of the operator lies in the closed left half-plane 
and the boundary part is diagonalizable (i.e. Jordan boxes corresponding to pure 
imaginary, including zero ones, eigenvalues are diagonal). The norm, with respect 
to which the operator is stable dissipative, exists if and only if the spectrum of the 
operator lies in the open left half-plane. 

Several stable dissipative (in their own norms) operators not necessarily are si- 
multaneously dissipative. To demonstrate that, consider operators represented by 



the matrices 



Ai 




A, 




Each of them is stable dissipative in its norm (due to the location of the spectrum) 

But 

/ -2 3 

Ai + A2 = 

V 3 -2 

The spectrum of the operator {Ai + A2) contains the point A = 1 which does not 
belong to the closed left half-plane. Thus, the operator {A1 + A2) is not dissipative in 
any norm. By lemma 6 the operators Ai and A2 are not simultaneously dissipative. 

The problem to find out necessary and sufficient conditions of simultaneous dis- 
sipativity of an arbitrary (even finite) family of operators seems to be wery difficult. 
Nevertheless, one can obtain some sufficient conditions imposing different constraints 
on the operators. Obtain the sufficient condition of simultaneous dissipativity of the 
family generating a solvable Lee algebra. Remind [13] that a family of matrices gen- 
erates solvable Lee algebra if and only if all elements of this family are simultaneously 
reducible to triangular form (generally speaking in complex basis). 

Theorem 6. Let the family {Aa} be compact and generate solvable Lee algebra, 
and the spectrum of each operator A^ lies in the open left half-plane. Then {Aa} 
is simultaneously stable dissipative. 

Proof. First consider the case of complex space C". Consider matrices of the 
operators A^ in the basis where they are of triangular form. 

Let each matrix A^ have the form 



Mi' 



A. 



(a) 











,,(a) \{a!) 







(a) 



\ 







,,W .,W ., (^\ .,W \{a) 



Show the existence of such a set of positive numbers {ck}k=i that all Aa are stable 
dissipative in the norm 

\Zk\ 



max 

l<fe<n Cfc 



(16) 



(here Zk is the k-th coordinate of the vector z in the given basis), whose unit ball is 
the polycylinder 

\zk\ <Ck{k=l,...,n). (17) 

If {ek}k=i is the considered basis, then, evidently, norm (16) coincide with the 
/°^-norm with respect to the basis {cfc/efc}"j=i in the norm (16): 

Re aii + Y,—\aij\ < {i = 1, . . . ,n). (18) 



j¥=i 



. c 



For the matrices A^ the conditions (18) look hke this: 
Re aS") 

C2 



Ke A^ ^ H |/i„i I 



,(°)i 



C?i— 1 I (a) I 



< 0; 

< 0; 

< 0. 



(19) 



Suppose /i = sup„;j_^; Ia^m h ^ = "Supo,^^, Re A^ . From the conditions of the 
theorem follows that < A < +oo, < /x < +oo. To fulfil (19) for all Aa, it is 
sufficient that the inequalities 

(ci + . . . + Ck^i)i2 < CfcA (A; = 1, . . . , n); ci > (20) 

be satisfied. 

Show the solvability of system (20). Let ci = 1. Choose the others Ck so that 

C2 > yu/A; C3 > (1 + C2)/i/A; . . . ; 

C„ > (1 + C2 + . . . + C„_i)/i/A. 

Then the inequalities (20) are satisfied, i.e. all operators Aa are stable dissipative 
in the norm (14). 

Let now operators A^ act in the space M". In usual way complexify M" and 
the family Aa- Then, as it has been described above, construct a cylinder (17). 
Intersection of (17) with the initial space M" produce a ball of the norm in which 
all Aa are stable dissipative. The theorem is proved. 

If instead of stable dissipative operators one considers dissipative operators, then 
the analog of theorem 6 is not true, starting from real dimension 4. Let 

V 2i y V i 

Each of the operators Ai^2 is dissipative in its norm. The finite family is compact, 
the matrices Ai and A2 generate solvable Lee algebra. Nevertheless 



A1 + A2 



The only eigenvalue of the operator {Ai + A2) is pure imaginary, with the matrix 
of this operator representing (up to a constant factor) non-trivial Jordan box. That 
means it is not dissipative in any norm, i.e. Ai and A2 are not simultaneously 
dissipative. To obtain a real example, one has to make the matrices Ai and A2 real: 

/O-ll 0\ /0-2 1 0\ 




A- 



R 



10 1 







Ai 











To keep true the statement about simultaneous dissipativity for nonstable dissi- 
pative operators, it is sufficient to strengthen the requirement of solvability up to 
nilpotency. Remind [13] that for each linear operator A in the space E the operator 
ad A in L{E) is defined: 

{ad A)B = AB- BA. 

The family {Aa} generates the nilpotent Lee algebra if and only if there exists 
such a number m & N that for any set of {A^^}^^ (among the elements of wich 
there may be the same ones) and for all a: 

m 
nM^aJ^a = 0. (21) 

fc=l 

Nilpotent Lee algebra is always solvable. Commutative Lee algebra is nilpotent 
(for it m = 1) and solvable. 

Theorem 7. Let the family {A^} be finite and generate nilpotent Lee algebra, 
and for each operator A^ exist a norm with respect to which it is dissipative. Then 
{Ak} is simultaneously dissipative. 

Proof. Without loss of generality one can assume that among the operators Ak 
there are no scalar ones (if A = al, where Re a < 0, then A is dissipative in any 
norm) and exists at least one operator (denote it Ai), among eigenvalues of which 
there are pure imaginary (otherwise we are under the conditions of theorem 6). 

First assume that A^ operates in C". We prove the theorem by induction on 
dimension of space. In dimension 1 the statement of the theorem is trivial. Show 
that one can expand all the space C" into a direct sum of two non-trivial subspaces 
invariant with respect to all At- Since in both of them the conditions of the theorem 
(for corresponding restrictions of {^fc}) are satisfied, then to complete the proof one 
has to use lemma 6. 

Let A be an imaginary eigenvalue of Ai; E' be the corresponding to A eigen- 
subspace (by virtue of diagonalizability of boundary part of Ai it coincides with 
whole corresponding root subspace); E" be the sum of root subspaces corresponding 
to all the others eigenvalues of Ai. Evidently, C" = E'®E" (the sign © means direct 
sum); E' ^ C", otherwise the operator Ai is scalar. Show the invariance of E' and 
E" with respect to all A^. 

Let X G E' . Then 

Aix = Xx. 

On the other hand, in accordance with (21) there exists such m E N that 
{ad Ai)'^Ak = for all k and 

{Al - XirAkX = 0. 

A more general fact is true: if Ax = and {ad A)"^B = 0, then A™'Bx = 0. For 
771 = the fact is obvious. Let that be true for m = r. Assume 

Ax = 0; {ad Ay+^B = 0. 

Then {ad AY{ad A)B = 0, and according to the inductive hypothesis 
A''{ad A)Bx = 0. But A''+^Bx = A'{BAx + {ad A)Bx), i.e. A'+^Bx = 0, as 
was to be proved. 

As a consequence of coincidence of E' with the whole root subspace, correspond- 



i.e. AiAkX = XAkX, AkX G E'. 

Show now the invariance of E". Let {ej}"'^i be the Jordan basis of the operator 

Ai with E' being corresponded to the vectors {ejYjLj^. One has to show that for 
any j less then ji or more then j2 the coordinates of A^Cj with the numbers from ji 
to J2 with respect to the assigned basis are equal to zero. Let it be not so and exist 
such j' that Cj' G E", but the ji-th coordinate (ji < jo < J2) of the vector A^ej' is 
a 7^ 0. Write it like this: 

AkCj' ^ . . . ^ oCjQ. 

Let Cj/ be an eigenvector of Ai corresponding to the eigenvalue fi ^ X. Then 

{ad Ai)Akej' = Ai{. . . + acj,,) — /i(. . . + acj;,) = . . . + (A — fi)aejg. 

Verify that 

{ad Ai)"'Akej' = . . . + (A - fi)'^aej^y 

For 771 = it is obvious. Let it be satisfied for m = r. Then 

{ad AiY^^AkCy = {ad Ai){ad AiJ AkCji = 

= Ai{ad AiYAkCf - {ad AiY AkAiCf = 
= Ai{. . . + {\ — nYa'^ja) '~ fJ'{(id AiY Ak&j' = 

= . . . + (A — ^lY^^^ja ~ /^(-^ ~ fiYdGjo = 

= . . . + (A - fiY^^o.^jo^ 

i.e. that is true also for m = r + 1, and, hence, for all m & N. 
Thus, 

{ad Ai^Akey = . . . + (A - /i)'"aej„ ^ 

for any m E N, which contradicts (21). 

Let now Cj' be a root (but not eigen) vector, corresponding to the eigenvalue /x, 
with the jo-th coordinate of the vector AkCj'-i equal to 0. Then 

{ad Ai)Akeji = Ai{. . . + acj^) — Ak{eji_i + jj^eji) = . . . + (A — fi)aejQ. 

Analogously 

{ad AiY'Akaj> = . . . + (A - /x)™aej„ ^ 

for any m E N , which contradicts (21). 

Since the sequence of basis vectors belonging to the root subspace begins with 
the eigenvector, the required statement for complex space is proved. 

The transfer onto the case of real space can be done in the same way as in the 
proof of theorem 6 (the ball of corresponding norm in the copmlexified M" intersects 
with R"). The theorem is proved. 

From theorems 6 and 7 follows, in particular, that a finite (compact) commutative 
family consisting of operators dissipative (stable dissipative) in their own norms is 
simultaneously dissipative (simultaneously stable dissipative). 



2.2 The Mass Action Law and Dissipative Mechanisms 

Some constructive conditions of simultaneous dissipativity can be obtained for finite 
families of operators of rank 1. The problem of the absence of Moore effect in the 
system constructed in accordance with the Mass Action Law (MAL) is reduced to 
the problem on simultaneous dissipativity of such operators. 

MAL systems appear from mathematical description of systems of chemical and 
biological kinetics and in some other problems. To the considered process is assigned 
an algebraic object, called reaction mechanism and having the form: 

arlAi + . . . + arnAn -^ /Srl^l + • • • + PrnAn, {r = l,...,d). (22) 

Speaking in terms of chemical kinetics, the reaction mechanism is a list of stoi- 
chiometric equations of elementary reactions (22). In this case Ai, . . . ,An are the 
substances taking part in the reaction; ari, (3ri are the non-negative integers called 
stoichiometric coefficients and showing in what amount the particles of Ai enter into 
the r-th elementary reaction as the initial substance (a^) or product (/9rj)- The fol- 
lowing notations are accepted: •jri = Pri — c^Hi Ir is the vector with the components 
7rj (i = 1, . . . , n) - so-called stoichiometric vector of the r-th elementary reaction. 

In accordance with MAL [14,15], to the mechanism (22) corresponds the following 
system of ordinary differential equations: 

-Tr = 2^ IriWr (23) 

at j.^^ 
where Cj(t) is the concentration of substance Ai at the moment of time t > 0, 

n 
Wr = K (t) Y\_ C^f^ 

i=i 

is the rate of the r-th elementary reaction, continuously depending on time. In 
particular, if reaction proceeds under constant external conditions, then kr{t) = 
const (r = 1, . . . ,d) and kr is called rate constant of the r-th elementary reaction. 
In the latter case (23) represents an autonomous system with polynomial right sides. 
Let L be a linear envelope of the family {'~fr}f=i- li L ^ R", then there exist such 
Oj (i = 1, . . . , n), not all equal to sero, that for all r = 1, . . . , d the equalities 

n 

Y^ ai'Jri = 
i=l 

are satisfied, from which for system (23) follows that 

n 

Y^aiCiit) = const (24) 

Relationships (24) are called stoichiometric conservation laws. If all aj are posi- 
tive, then the corresponding stoichiometric law is called the positive conservation law 
[15]. In MAL positive conservation laws takes place rather often (but not always). 

As it is known [15], balance polyhedrons are intersections of affine subspaces of 



to (23) convex sets (one can find the proof of tlieir positive invariance in [15]). If 
there exists at least one of positive conservation law, they are compact. 

The question arises: under what conditions does the norm exist in R" according 
to which the system (23) is constricting in all balance polyhedrons and independent 
of rate constants? 

Definition 8. Mechanism (22) is called dissipative, if for system (23) there 
exists a norm, constricting in all balance polyhedrons irrespective of rate constants 
(in other words, the constricting norm depends on the mechanism only). 

We use the notation M^i for the operator in M", represented by the matrix, in 
the i-th column of which there are components of the vector 7^, and on other places 
- zeros. The subspace L is invariant with respect to all M^ [15]. The notation M^j 
stays for restriction of Mri on L. 

Theorem 8. Let for mechanism (22) exist at least one positive conservation 
low. This mechanism is dissipative if and only if the family {M^j : a^ > 0} is 
simultaneously dissipative. 

Proof. Sufficiency. It is known [15] that the Jacobi matrix J^ of system (22) at 
the point c, whose coordinates are positive, has the form 



Ci 



Jc= J2 ar^^Mri. (25) 

ari>0 



Matrices Jc belong to the convex cone produced by the family {Mri\ a^ > 0}. 
Besides, the difference of any two solutions (23) from one balance polyhedron belongs 
to the subspace L. Under the conditions of lemma 3 and theorem 2 obtain the 
existence of constricting norm in the subspace L. It can be expanded onto all M". 

Necessity. Matrices Mri {c(ri > 0) belong to the closure of the family of matrices 
Jc for arbitrary non- negative vectors c and rate constants kr. To prove this, first 
let consider the case when Cj {j = 1, . . . ,n) and kr are fixed and all ki{l ^ r) tend 
to zero. In the limit in (25) only the sum for given r is left. Further on, fix all 
Cj > (j 7^ i) and let q tend to zero, changing kr so that the equality ariWr/ci = 1 
holds true. Then all the terms except one tend to zero and in the limit we obtain 

Mri. 

Thus, the matrices M^ (a^ > 0) belong to the closure of the family of restrictions 
of the matrices Jc on the subspace L. Hence, according to lemma 3 and theorem 2 
the necessity follows. The theorem is proved. 

Note that matrices Mri represent matrix-columns (in each matrix there is only 
one non-zero column) and that means that the rank of each of them is equal to 
unity. We come to the problem of simultaneous dissipativity of the finite family of 
operators of rank 1. 

Note that dissipative mechanisms of reactions were studied in details in [12]. In 
particular, some classes of dissipative mechanisms are pointed out and all dissipative 
mechanisms for n = 3, I]f=i ciri < 3, J2^^i /5ri < 3 (r = 1, . . . , d), C1 + C2 + C3 = const 
enumerated. 

In the next subsection are obtained necessary and sufficient conditions of simul- 
taneous dissipativity of the operators of rank 1 in R^ (corresponding to the case 
dimL = 2) and some sufficient conditions of simultaneous dissipativity of matrix- 
columns. 



2.3 Constructive Conditions of Simultaneous Dissipativity of 
One-Dimensional Operators 

Before consideration of simultaneous dissipativity of operators of rank 1, find out 
what can be said about dissipativity of one such operator. From necessary and 
sufficient conditions (see the paragraph after lemma 6) follows that the norm in 
which the given operator of rank 1 is dissipative exists if and only if it has a negative 
eigenvalue. 

Positive semi-trajectories of system (3) corresponding to the initial condition 
a;(0) = xq are in this case rectilinear segments parallel to the image of A and 
connecting xq with Ker A. Operator A of rank 1 is dissipative in the given norm if 
and only if for any point x {\\x\\ = 1) there exists such e > that ||a; + eAx\\ < 1. 
It means that the negative number belongs to the spectrum of A, and the image of 
A is orthogonal to its kernel (in the given norm, the subspace E2 is orthogonal to 
El if \\x + y\\ > \\x\\ for any x G Ei,y G E2 [7]). 

Let now be given a family {M^}]^^^ of operators of rank 1 in M". Each of them 
can be represented in the form (■ ;-?/'fc)y9fc, i.e. M^x = {x ;ipk)fk where (■ ; ■) is 
the standard scalar product in M". The vectors (fk and ipk are determined by the 
operator Mk unambiguously (up to scalar factors). Let Xk = {ipk','ipk), i-e. A^ is an 
eigenvalue of Mk (either it is the only non-zero eigenvalue, or 0, if the operator Mk 
is nilpotent). As it has already been mentioned, for simultaneous dissipativity of 
{Mk} the conditions 

Xk<0 {k = l,...,m) (26) 

are necessary. 

Assign to each operator Mk the projector Pk projecting parallel to the image of 
Mk on the kernel of Mk- It is easy to see that Pk = I — MkjXk- By virtue of the 
above-mentioned condition of dissipativity of the operator of rank 1 in the given 
norm the operator Mk is dissipative in some norm if and only if Pk is constriction 
in this norm. 

All Pk can be constrictions in one norm if and only if all products of the form 
11^=1 Pkj {q & N is arbitrary; kj G {1, . . . , m} and they are not necessarily different) 
are jointly bounded. We come to the following conclusion. 

Lemma 7. The family {Mfc}^^ of operators of rank 1 is simultaneously dissi- 
pative if and only if the conditions (26) are satisfied and all products of the form 
11^=1 Pkj (g G iV is arbitrary; kj G {!,..., m} and they are not necessarily different) 
are jointly bounded. As a constricting norm one can take 

o>' ^27) 

=N,1< kj<m " j^i ^ 

where ||.||o is any norm in M". 

Proof. All statements of the lemma, exept the latter, follow immediatily from 
the above reasonings. Further on, if all products n^=i Pk are jointly bounded, then 



sup {|a:|o, 


{[[Pk,)x 


qeN,l< kj<m 


.7 = 1 



sup {||a;||o, Ikn^fej)^ nJ' < °° 

q£N,l< kj<m j—^ 



0' 



for each x G M". This expression possesses all properties of norm and all operators 



From lemma 7 follows a simple consequence. 

Corollary 1. If all ipk are collinear (images of Mk coincide) or all ipk are collinear 
(kernels of M^ coincide) and {ipk, "ipk) < for all /c = 1, . . . ,m, then the operators 
Mk {k = l,...,m) are simultaneously dissipative. As corresponding constricting 
norm one can take 

sup {||a;||o, ||Pfca:||o}- 

qeN,l<kj<m 

To demonstrate this, it is sufficient to note that in these cases 

Q 

n Pi^j = ^fci 

or 

q 

Y[Pk, = Pk,, 

i=i 
respectively. 

Remark 10. If not all (pk are collinear, then as a norm in lemma 7 one can take 

sup \\{Y[Pk,)x . (28) 

q&N,l<kj<m" j—i " 

The criterion established in lemma 7 is not constructive. Constructive criteria 
of simultaneous dissipativity of finite family of operators of rank 1 in R" have been 
obtained only at n = 2 (for arbitrary n there exist sufficient conditions for one class 
of operators; they are given at the end of the section). Pass to the consideration of 
the case n = 2. 

Consider the family {Mk}^^i of the operators of rank 1 in M^. As before, represent 
each operator Mk in the form (■ ; ipk)fk- Let first m = 2. 

Lemma 8. The operators Mi = (■ ; ipi)^pi and M2 = (■ ; '?/'2)v^2 are simultaneously 
dissipative in M^ if and only if the condition 



(<^i;^2) ■ (^2;^i) 



< 1 (29) 



(v?i;^i) ■ (^2; ^2) 
is satisfied together with the conditions 

i^uiJi) <0; (^2;V^2) <o. 

As a corresponding constricting norm one can take 

||a;|| = max{||a;||o, ||Pia;||o, ||P2a;||o, ||PiP2a:||o, ||^2^ia;||o}- (30) 

Proof. In R^ the projectors Pi and P2 have rank 1 and are represented in the 
form 

Pi = (■ ;^i)xi; P2 = (■ ;^2)x2, 

where r/i, 772, Xi? X2 are some vectors in R^. 

The operators nj=i Pkj are bounded when the spectrum of the operator (-P1P2) 
lies on the segment [—1; 1]: 



In a standard orthonormalized basis Pk acts like this: 

PkX 




where 

denotes the vector with the coordinates a*^^-* and a^'^\ Hence 

(v^i;V'2) 



(^i; X2) 
(xi;^2) 



(v^2; V^2) 
(^2;^i) 



(v'l;^!) 

i.e. condition (31) takes the form (29). 

To complete the proof, use lemma 7. To check a possibility of choosing corre- 
sponding norm in the form (30), note that 

iPiP2yPi = ivi;x2r-iv2;xir-Pi; 
iP2PiYP2 = ivi;x2Y-iv2;xiY-P2 

for any r & N. It means that with the account of (31), in (27) one can restrict 
oneself to finite number of products. The lemma is proved. 

Remark 11. If v^i and ip2 are non-collinear, then as required norm we can take 

max{||Pia;||o, ||P2a;||o, \\PiP2x\\o, \\P2Pix\\o}. 

This follows from remark 10. Then the ball of the norm is determined by the 
inequalities 

'1 1 1 



{x; r]i)\ < min< 
{x; 772)1 < min< 



IXillo' l(Xi;^2)| ■ IIX2II0, 
1 1 



IX2II0 I(x2;m)| ■ llxilloj 

i.e. it is parallelogram. 

Also note that for simultaneous dissipativity of a family the dissipativity of each 
operator from convex envelope of the family is insufficient. To see this, consider the 
operators represented by the matrices 

/-I 1\ / 

Ml = ; M2 = 

Each of them is dissipative in its norm. It is easy to show that spectrum of any 
non-trivial convex combination of Mi and M2 lies in open left half-plane. Neverthe- 
less 

(yi;V^2) • {^2; A) ^ _2 



Reasoning like in proof of lemma 8, it is easy to obtain a criterion of a simulta- 
neous dissipativity for arbitrary m. The result is a set of conditions of the form 



ifhi^k) <0 {k 



,m 



(v'fcii^fe) • (v'fei^fcs) ■•••■ i^k,;ip, 



ki, 



<1, 



(32) 
(33) 

(•y^fei; V'fei) ■ (v'fea; V'fca) ■ • • • ■ (v^fc,; V'fcJ 

where {/cj-j^^^ is a set of different numbers from 1 to m, and inequalities (33) holds 
for all such sets. The number of conditions has the order 0{{in — 1)\) and for any 
large m testing of these conditions becomes unrealizable. It turns out, however, that 
among inequalities (33) there are dependent ones and the number of conditions can 
be reduced. 

Theorem 9. Let the vectors il)^ {k = 1, . . . ,m) lie in one half-plane clockwise. 
Then the family of operators {M^j^j^ where M^ = {■;ipk)Vk is simultaneously dis- 
sipative if and only if the vectors (fk {k = ^, ■ ■ ■ ,itl) lie in one half-plane clockwise 
and the conditions (32) and the followings ((34), (35)) are satisfied: 



{ipk;ipk+i) ■ {(pk+h-ipk) 



< 1 



(k 



{Vk;'ipk) ■ {(Pk+ui^k+i 

1, . . . , m with iprn+l = -V^i; ■ipm+1 

(<^i;^2) ■ (^2; ^3) ■•••■ (^m;V^i) 



(34) 



(v^i;V'i) ■ (v^2;^2) ■ ••• ■ ivm;i^m] 



< 1; 

< 1. 



(35) 



(V9i;V'l) ■ (9^2; ^2) ■ ••• ■ iVm;^m) 

The corresponding norm can be chosen polyhedral (a norm, whose ball is polygon). 
Proof. Necessity. Let Ik be the kernels of the operators M^ (i.e. straight lines 
orthogonal to ipk)- Straight lines Ik divide the plane into 2m sectors. If among 
the vectors ipk there are collinear, then some sectors are singular, but this does not 
change the further reasonings. In each sector G and for each p E {!,..., m} 

sign{xi]i)p) = sign{x2;^p) 

for all Xi G int G, X2 € int G. 

Let Gr be a sector lying between corresponding rays of straight lines /,. and /,.+i 
(where Im+i = ^i)- It is enough to consider the sectors {Gk}^=i into which one 
half-plane is divided, since for sectors lying in vertical angles to Gk the reasons are 
the same. 

Note that by inequality (32) for each operator Mk the projector Pk is determined, 
which operates in each sector Gk as a projector in the direction Vkr = sign{x; ipk) ■ 
(fk {x G Gr) onto the straight line Ik- 

A norm with respect to which all Mk are dissipative exists if and only if there 
exist a convex body Q symmetrical with respect to and positively invariant with 
respect to all systems of the following form 



dx 
Itt 



^hk{t){x]^k)Vk, 



(36) 



where hk{t) is any function piecewise continuous and non- negative for t > 0. The 
sufficiency is evident (suppose hk{t) = l,hj(t) = for j 7^ A; and come to dissipa- 
tivity of Mk with respect to Q). To prove the necessity, it is sufficient to make an 
estimation analogous to that made in the proof of theorem 1: 

J m 

-Mt)\\Q = NQ[x{t), J2 hk{t)Mkx{t)) < 
"■'^ fc=i 

m 

<lQ{T.hkit)Mkxit))-\\xit)\\Q<0. 

k=l 

Here ||.||q is a norm whose unit ball is Q. 

Since (x; ipk)'fk = \{x; ipk)\ ■ ^kr at a; G Gr, then (36) can be rewritten as follows: 

dx "" 

-^ = Y.ykitWr (37) 

where yk{t) is piecewise continuous and non- negative for t > 0. Thus, it is sufficient 
to construct such a polygon W that from each point of its boundary dW all the 
vectors Vkr sue not directed into the exterior of W. Then one can take 

Q = co{WU{-W)}. 

Let (37) have at least one unbounded solution, whose positive semi-trajectory 
lies inside one of sectors. Then (36) has an unbounded solution, i.e. the operators 
Mk are not simultaneously dissipative. 

The notation C{vkr} is used for a convex cone produced by {vkr}^=i- 

Let this cone concide with M^ at least in one sector Gr (i.e. the vectors generating 
it do not lie in one half-plane). Then as yk{t) one can choose such constants that 
''J = YJk=iyk'^kr G Gr, and then, drawing a ray from the point x^ G intGr in the 
direction of v, obtain a positive semi-trajectory of unbounded solution (47) lying 
inside Gr- 

Thus, for simultaneous dissipativity of {Mk} it is necessary to satisfy the condi- 
tions 

C{t;fcJ ^M^ (A;= l,...,m). (38) 

If C{vkr} in some sector is a half-plane, then it must contain the vertical angle to 
Gr - Gr (and thus intersect with Gr only at zero); otherwise (37) has an unbounded 
solution. For each sector Gr consider the boundary of the cone C{vkr}- It consists 
of two directions. Show that for Gj it is Vjj^vi^j^i-) j. It is sufficient to show that for 

J = 1- 

Let t;i^i and 1^2,1 be collinear and oppositely directed. Then to satisfy (38) it is 
necessary tha the other Vki lie on one side of the straight line, stretched on t;i^i. But 
if i^i^i and t;2,i are non-collinear, then all other Vki can be expanded in terms of the 
basis 1^1,1,1^2,1- 

Let, for example, t;3,i = cit;i,i -|- C2i'2,i, and 113,1 be collinear to one of the basis 
vectors (for example, iii,i; the case with i;2,i is considered analogously). Then C2 = 0. 
If Ci > then 173,1 G Cl^i, 1,112,1}. Let Ci < 0. Then to satisfy (38) in Gi it is 
necessary for i^ii and 113,1 to be boundary directions in C{i;fci}. Since i'fc,(/+i) = '^kh 



C{vki}, 1^2,2 = — '^-'2,1 £ C{vk2}- It means C{vki} and {vk2} represent half-plane 
whose join is all M^, what is impossible. That means ci > 0. 

Let now v^^i be non-collinear neither to i^i^i nor to t;2,i. If Ci < 0, C2 < 0, then in 
Gi (38) is not satisfied. If Ci < 0, C2 > 0, then in G2 there Vs^2 = CiVi^2 + (—£2)1^2,2, 
i.e again (38) is not satisfied. Analogous reasonings hold for the case ci > 0, C2 < 0, 
i.e. the only possible case is Ci > 0, C2 > and therefore v^^i G C{t;i^i,i;2,i} (where 
C{x,y} is a convex cone, stretched on the vectors x and y). 

The case is left when the directions Wi^i and t;2,i coincide. 

Without loss of generality one can assume non-collinearity oi v^^i and t;i^i. Then 
772,2 and 1^3,2 are boundary directions in C{vk2}- Consequently, t^i^i G C{—Vi^i,V3^i} 
i.e. the directions ^3,1 and 1^1,1 coincide contrarily to the assumption. It means that 
if Vrr and V(r+i),r are co-directed, all Vkr are collinear, i.e. all ipk are collinear. In 
this case the directions Vrr and V(^r+i),r are also boundary. 

We call the obtained fact the boundariness condition. 

Since all ipk lie clockwise in one half-plane, then it is easy to check that in sector 
Gm either all {x; ipk) > for all k or {x; ipk) < for all k. Thus, by virtue of 
(32), all ifk lie in one half-plane. From the boundariness condition follows that 
(fk G C{ipk-i,^k+i}) i-e vectors (pk are arranged either clockwise, or anti-clockwise. 

Let, for example, i;i,i = ipi (the case i;i,i = —ipi is considered analogously). 
Then i;2,i = ^2 lies in the half-plane bounded by the straight line stretched on vi 
and containing Si. Therefore the direction from ipi to ip2 in the half-plane containing 
all (pk is the same as from ipi to il)2-, i-e. clockwise. 

The necessity of the other conditions is obvious, since (34)-(35) is simply a part 
of conditions (33). 

Sufficiency. Let the family {ipk}^=i be arranged clockwise in one half-plane and 
the conditions (32) and (34)- (35) be satisfied. Assume that among ipk there are 
non-collinear vectors, and among ipk there are no collinear ones. 

The condition of clockwise arrangement of (pk in one half-plane means that 
the angle (counted from (pi clockwise) between (pi and the vectors v^i, v^2, • • • ,V^m, 
^m+i = —^i monotonously increases from to vr. Taking into account that the an- 
gle between (pki and {—(pk2) is the angle between (pki and ipk2-i taken with opposite 
sign, it is easy to conclude that systems {i^fcr} (in each sector) lie in one half-plane 
and are arranged clockwise (to avoid exiting from corresponding half-plane we start 
counting in sector Gr from W(r+i),r)- 

From conditions (34) follows that in each sector there is a " convex configuration" , 
i.e. there is vector x G Gr, representable in the form 

m 

•^ / J CkVkri 

fc=l 

where all Cfc > 0. 

It means that if from one point x G int Gr one draws segments a and h in the 
directions of Vrr and V(^r+i),r up to the crossing with Ir and Z^+i, respectively, then 
these segments together with the segments connecting with the point of crossing 
a with Ir and h with Z^+i, respectively, form a convex polygon (if Vrr and V(j.+i),r are 
oppositely directed, it will be a triangle, and if they are non-collinear - a quadrangle; 
as we have seen before they cannot be co-directed). 

Due to the same orientation of {pk} and {ipk} all the other Vkr are directed (from 



Fix now the point Xq G h {xq 7^ 0) on the boundary ray of sector Gi (actually, 
one can begin from any straight line Ik', we begin from /i). Due to the boundariness 
condition either direction from xq on I2 goes into sector Gi, or direction from xq on 
Im goes into G„i- 

If one and only one of these statements is true, continue moving in the correspond- 
ing direction (to the neighboring straight line) till the direction on the neighboring 
straight line goes into the neighboring sector. In other words, move from /^ to Z^+i 
in the direction parallel to ipr, if this direction goes into sector Gr (or, into Gr-i, 
respectively). As a polygon W mentioned after (37) one should take a polygon 
formed by the segments which we moved along, and the segments of those straight 
lines on which the movement broke (if exit on the initial ray did not occur, in our 
case it is a part of /i corresponding to Gi, then it is a segment connecting xq with 
0, and a segment of that straight line on which the movement broke, connecting the 
point of breaking with zero; if exit on the initial ray occured, then it is a segment 
connecting xq with the point of exit). 

If both statements are satisfied, then as W one can take a join of two such 
polygons formed in moving to both sides from Xq. 

This algorithm is easy to check proceeding from boundariness conditions, "convex 
configuration" , and (35) (the latter condition means that if exit on the initial ray 
occured in moving in either side, then the point of exit is no farther from the 
beginning of coordinates than the initial point; in particular, if the point of exit 
coincides with the initial point, then the formed polygon can be taken as W). The 
ball of the sought for norm is a polygon. 

If some of ipk are collinear, then some sectors Gk are singular. This, however, 
does not change the results. The reasonings are analogous to the case when among 
ipk there are no collinear vectors. The only difference here is the following: some 
straight lines Ik correspond to several directions {v^jj^i;,,,. Then in constructing W 
one needs to move along ipko- 

In the case when all ipk (or all ipk) are collinear (see corollary 1), all the same one 
can regard that {^pk} and {ipk} have the same orientation, starting from (32). 

Conditions (34)-(35) are satisfied in this case. The norm can be chosen polyhe- 
dral, if one chooses a polyhedral norm as ||.||o in (30). The theorem is proved. 

Remark 12. One can obtain the arrangement of vectors ipk required by the 
conditions of theorem 9 by renumbering vectors and (if it is necesary) changing 
signs of some of them. 

Thus, the problem of simultaneous dissipativity of a family of operators of rank 1 
in M^ is solved completely. The number of conditions to be checked now, in contrast 
to (33), is only of the order 0{m). 

With theorem 9 one can study the MAL mechanism on dissipativity (and, re- 
spectively, on the absence of IME). For example, let the mechanism be 

Ai -^ A2, Ai -^ A3, A2 -^ Ai, A2 -^ A3, 
3^2 ^ Ai + 2A3, 2Ai ^ A2 + A3, 2^2 ^ Ai + A3, 
2A3 ^ Ai + A2, 3Ai ^ A2 + 2A3, 3A2 ^ 2Ai + A3, 

Ai + A2 ^ 2A3. (39) 

This mechanism possesses positive conservation law ci + C2 + C3 = const. The 



Obviously, dimL = 2, and one can use theorem 9. Writing matrices M^^ and using 
theorem 9, let make sure that mechanism (39) is dissipative. The corresponding 
norm in the subspace L has the form 



l|c|| 



|Cl| 



|C2|- 



It can be expanded onto all M^, for example, in this way: 

l|c|| = |ci| + |C2| + |Ci + C2 + Cs\. 

To complete the section, consider the question of simultaneous dissipativity of the 
finite family of operators of rank 1 of special form in M" for arbitrary n. Namely, we 
consider operators represented by matrix-columns. Let obtain sufficient conditions 
of simultaneous dissipativity of such operators. 

Let the basis {ckj^^i and the norm 



\x\ 



'^Pk\Xk\ 
k=l 



(40) 



be given in M", where pk > {k = 1, . . . ,n), Xk is the k-th coordinate of vector x in 
the basis {cfc}. Norm (40) coincides with /^ norm with recpect to the basis {ck/pk}- 
Therefore, the necessary and sufficient dissipativity conditions of the operator A 
represented by the matrix (ay)". ^ according to remark 9 have the form 



Piaa + "^Pjlajil < (i = 1, 



,n). 



(41) 



Let now there be a family of operators, represented by the matrix-columns Aki^ 
{k = 1, . . . ,n; Ik = 0, . . . ,rfc), where Aki,, is the Ik-th matrix with non-zero A;-th 
column: 

/ ... afl:^ ... \ 



AkU 







4fc^ 







Mk) 



(42) 



Vo ... cC' ... oy 

Coming from (41), write dissipativity conditions of all operators in norm (40) 
with some constants pk'- 



PfcoL +EpjI 



^jk I 



< (/c = l,...,n; lk = 0,...,rfc). 



(43) 



jVfc 



Theorem 10. If the system of linear inequalities (43) complemented by the 
inequalities 

Pfc>0(A; = l,...,n) (44) 

has a solution, then the family of operators represented by matrices (42) is simulta- 
neously dissipative. 

Proof. Solvability of the systems (43)- (44) means the existence of positive con- 
stants Pk {k = 1, . . . ,n) ioi which inequalities (43) are satisfied, and that is dissipa- 



Thus, for simultaneous dissipativity of finite family of operators represented by 
matrices-columns the solvability of above written finite system of linear inequali- 
ties proves to be sufficient. To check solvability, one can use algorithms of linear 
programming [16]. 

Remark 13. The solution of the system (43)-(44) exists if there exists solution 
of the system of (n — d) linear inequalities complemented by inequalities (44) (where 
d is the number of those k for which r^ = 0; evidently < d < n — 1). To prove 
this, 



assume 



Ofcfe 



max a^^ , 



Ijk 



max 



Aik)\ 
^jk I 



(j 7^ k) 



{k = l,...,n). 



Consider the system 






n 



(45) 



Obviously, if the set {pk} satisfies the system (44)- (45), then it satisfies the system 
(43)- (44) as well. Numbers k for which r^ = are excluded. Therefore, in system 
(45) there are {n — d) inequalities. 

Remark 14. For n = 2 theorem 10 provides necessary and sufficient conditions 
of simultaneous dissipativity. To demonstrate that, note that for operator M^ of the 
considered form the vector ip^ (see the notation at the beginning of the subsection) 
is directed along one of the coordinate axes. Therefore (see the proof of sufficiency 
in theorem 9), if the family is simultaneously dissipative, then one can choose paral- 
lelogram as a ball of the corresponding norm, with vertices on coordinate axes, i.e. 
the norm is of the form (40). In the case of arbitrary n the conditions of theorem 
10 are already not necessary. To see this, let 



/ -1 \ 



Ai 



/ -1 \ 



A, 



0-10 

Vo 1 oy 



-10 

V 1 oy 

The system of linear inequalities 

-P1+P2+P3 < 0; 

P1-P2+P3 < 

has no positive solutions. Nevertheless, simultaneous dissipativity exists, since each 
of the operators is dissipative in its norm and ipi = (f2 (see corollary 1). 



Conclusion 

Let us resume. The infinitesimal Moore effect (IME) in the interval space for 
smooth autonomous system on positively invariant convex compact is studied. The 
local conditions of absence of IME in terms of Jacobi matrices field of the system are 



of the Jacobi matrices is established, and some sufficient conditions of simultaneous 
dissipativity are obtained. 

On the basis of the conducted analysis the reason of weak efficiency of interval 
stepwise methods is pointed out. The main reason is that to solve the problem of 
absence of IME in the system and to construct corresponding interval space one 
needs analysis of simultaneous dissipativity of Jacobi matrices of system and con- 
structing a constricting norm. The latter questions are rarely solved constructively. 
Besides, in sufficiently rich interval spaces (for example, in using standard intervals 
- rectangular parallelepipeds) IME is almost always present. One should, however, 
remember that the notion of the Moore effect in the work is treated sufficiently 
strongly. The final conclusion on the efficiency of stepwise interval methods can be 
drawn only after studying asymptotic Moore effect (AME). It should also be noted 
that there may be definitions of interval spaces, different from definition 1. 

Some particular classes of systems without IME and corresponding interval spaces 
are pointed out. These results can be used in solving by interval methods particular 
systems from the pointed out classes. 



References 

1. Moore R.E. Interval analysis- N.-Y.: Prentice-Hall, 1966. 

2. Kalmykov S.A., Sliokin Yu.L, Yuldasliev Z.Kli. Methods of interval analysis- 
Novosibirsk: Nauka, 1986.- 221 p. 

3. Chernousko F.L. Optimal guaranteed estimations of uncertainty by means of 
ellipsoides // Izv. AN SSSR.- Tekhnich. kibernetika.- 1980.- N 5.- P.5-10. 

4. Kracht M., Schoeder G. Zur Intervallrechnung in linear Raumen.- Computing- 
1973.- V.ll.- P. 73-79. 

5. Ratscliek H. Nichtnumerische Aspecte der Intervallarithmetik / / Interval Ma- 
thematics.- Berlin-Heidelberg: Springer- Verl., 1975.- P.48-74. 

6. Kuratovskii K. Topology. V.l.-Moscow: Mir, 1966.- 594 p. 

7. Belitskii G.R., Lyubich Yu.L Matrix norms and their applications- Kiev: 
Naukova dumka, 1984.- 151 p. 

8. Kantorovicli L.V., Akilov G.P. Functional analysis- Moscow: Nauka, 1977.- 
706 p. 

9. Lozinsky S.M. Error estimation of numerical integration of ordinary differential 
equations / / Izv. vuzov. Ser. mat.- 1958.- N.5.- P.52-90. 

10. Bylov B.F., Vinograd P.A., Grobman D.M., Nemytskii V.V. Lyapunov index 
theory and its applications to the problems of stability- Moscow: Nauka, 1966.- 
576 p. 

11. Verbitskii V.I., Gorban A.N. Simultaneously dissipative operators and their 
applications in dynamical systems- Krasnoyarsk, 1987.- 32 p. (Preprint /AS USSR, 
SB, Computing Center). 

12. Verbitskii V.I., Gorban A.N. Thermodynamical restrictions and quasi-termo- 
dynamicity conditions in chemical kinetics / / Mathematical problems of chemical 
kinetics / Ed. by K.I.Zamaraev and G.S.Yablonskii.- Novosibirsk: Nauka, 1989.- 
P.42-83. 

13. Burbaki N. Lie groups and algebras- Moscow: Mir, 1976.- 496 p. 

14. Volpert A.I., Khudyaev S.I. Analysis in the classes of discontinuous functions 
and equations of mathematical physics- Moscow: Nauka, 1975.- 394 p. 

15. Gorban A.N., Bykov V.I., Yablonskii G.S. Essays on chemical relaxation- 
Novosibirsk: Nauka, 1986.- 300 p. 

16. Golshtein Ye.G., Yudin D.B. New tendencies in linear programming .- Mos- 
cow: Sov. radio, 1966.- 524 p. 

17. Verbitskii V.I., Gorban A.N., Utjubaev G.Sh., Shokin Yu.L Moore effect in 
interval spaces // Dokl. AN SSSR.- 1989. V. 304, N 1.- P. 17-21. 

18. Bykov V.I., Verbitskii V.I., Gorban A.N. On one estimation of solution of 
Cauchy problem with uncertainty in initial data and rigt part // Izv. vuzov, Ser. 
mat.- 1991. N. 12.- P.5-8. 

19. Verbitskii V.I., Gorban A.N. Simultaneously dissipative operators and their 
applications // Sib. mat. jurnal.- 1992. V.33, N 1.- P.26-31. 

20. Verbitskii V.I., Gorban A.N. Simultaneously dissipative operators and quasi- 
thermodynamicity of the chemical reactions systems // Advances in Modelling and 
Simulation, Tassin (France): AMSE Press. 1991.- V.26, N 1.- P.13-21. 

21. Verbitskii V.I., Gorban A.N. On one approach to the analysis of stability of 
nonlinear systems and differential inclusions // Advances in Modelling and Analysis, 



22. Verbitskii V.I., Gorban A.N. Stability analysis and solution evaluation for 
nonlinear systems by "Jacobian fields" and Liapunov norms / / AMSE Transactions, 
Scientific Siberian, A, V. 4. Dynamics. Tassin (France): AMSE Press. 1992.- P.104- 
133. 



